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We employ the S — t- vr form factors obtained from QCD light-cone sum 
rules and calculate the B — )• -kIv^ width {i = e,//) in units of 
integrated over the region of accessible momentum transfers, < < 
12.0 GeV^. Using the most recent BABAR-collaboration measurements we 
extract jKbl = (3.50l[jJ^|^^ ±0.111^^.^) x IQ-^. The sum rule results for 
the form factors, taken as an input for a z-series parameterization, yield the 
g^-shape in the whole semileptonic region of i? — >■ ttIvi. We also present the 
charged lepton energy spectrum in this decay. Furthermore, the current sit- 
uation with B — TVr is discussed from the QCD point of view. We suggest 
to use the ratio of the B — t- vrrz^,- and B — )• ttIv^ {i = /i, e) widths as an 
additional test of Standard Model. The sensitivity of this observable to new 
physics is illustrated by including a charged Higgs-boson contribution in the 
semileptonic decay amplitude. 
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1 Introduction 



Currently, there is a tension between the two values of \Vub\ extracted from inclusive and 
exclusive semileptonic -B-decays involving b ^ u transition. While the inclusive analyses 
typically yield a central value of \ Vuh\ larger than 4 x 10~^, the exclusive determinations 
produce central values well below this. This tension is not significant; it ranges at a 
level of 3cr, but it already has created a significant amount of speculations concerning 
possible new physics effects. This is in contrast to the situation with \Vch\-, where both 
the inclusive as well as the exclusive determinations yield consistent values with an 
uncertainty of roughly 2% (for a review on |V"cb| and \Vub\ see [1]). 

The theoretical description of inclusive semileptonic B decays relies on the heavy- 
quark expansion which has reached a mature state. Still, the situation with 6 — t- -u 
inclusive decays is more complicated than with the dominant 6 — >■ c ones. In order to 
suppress the charm background in the inclusive 6 — )• u decays, severe phase space cuts are 
necessary, for which most theoretical methods cannot rely on the heavy quark expansion 
based on the local operator-product expansion. The heavy quark expansion for this case 
uses non-local matrix elements corresponding to the light-cone distribution functions of 
the B meson, which also appear e.g., in the radiative b ^ s decays. Due to this more 
complicated structure of the expansion, it is quite hard to estimate subleading terms 
in the heavy-quark expansion for 6 — )• u decays. In particular the inclusive method has 
been scrutinized for a missing systematic effect like e.g., "weak annihilation", however, 
no missing pieces could be identified yet. Still, it is believed that this method allows a 
determination of \Vub\ at a precision of a little better than 10%. 

As far as exclusive decays are concerned, heavy-quark symmetries restrict the form 
factors for the heavy-light 6 — t- u transitions much weaker than the ones for the 6 — )• c 
transitions. Hence, the determination of \Vub\ from exclusive decays such as -B — t- ti^D 
and B — )• plv requires a QCD calculation of the relevant hadronic form factors. State- 
of-the-art calculations do not rely on quark models any more, since the latter cannot be 
directly related to QCD. Already for many years, the form factors are obtained, on one 
hand, from lattice simulations and, on the other hand, from QCD sum rules. The two 
approaches are complementary: While the lattice techniques can calculate close to the 
maximal leptonic momentum transfer g^, the QCD sum rule approach works best for 
small q^. Extrapolating both predictions to the full phase space yields consistent results 
and hence there is some confidence that the form factors, in particular, for B — >■ ttIv are 
known with an uncertainty of 10-15%. 

Currently, B — )• 7r/f£ is the most reliable exclusive channel to extract \Vub\- There is a 
steady progress in measuring the branching fraction and ^^-distribution for I = fi,e (see 
[21 El H] for the latest results). The hadronic vector form factor f^^{q^) and its scalar 
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counterpart f%^{q^) relevant for this decay are defined as0 



n2 ^2 



+ r^M') - g^, (1) 

where /gj^CO) ~ /sttC^)- "^^^ most recent lattice QCD computations with three dynam- 
ical flavours [5l |6] predict these form factors at > 16 GeV^, in the upper part of 
the semileptonic region 

< g2 < (tub - m^f ~ 26.4 GeV^, with an accuracy reaching 
10%. There are also recent results available [7J in the quenched approximation on a fine 
lattice. QCD light-cone sum rules (LCSR) with pion distribution amplitudes (DA's) 
allow one to calculate the B ^ n form factors [SI [9l [lOl [TTl [12] at small and intermediate 
momentum transfers, < < q^ax^ where the choice of q^ax varies between 12 and 16 
GeV2. 

The main goal of this paper is to present an updated LCSR prediction for the width 
of S — )• 7r^f£ (i = e, ii) which is then used to extract \Vub\- This work complements |I2] 
Bn and /^^ 



where LCSR for ft and /R were rederived, employing the MS scheme for the virtual b 



quark in the correlation function. In [12], the shape of the form factor /^^(g^) predicted 
from LCSR was fitted to the earlier BABAR measurement [13] of the ^^-distribution in 
B — nlvi. In this way, some input parameters of LCSR were constrained, allowing one 
to decrease the theoretical uncertainty of the value /^^(O), the main prediction of [T2] . 

In this paper, we follow a different strategy. The intervals of the Gegenbauer moments 
of the pion twist-2 DA are constrained using the LCSR for the pion electromagnetic 
(e.m.) form factor at spacelike momentum transfers. The calculated form factor is then 
fitted to the available experimental data on this form factor. We also slightly update 
the other input parameters, and recalculate the form factors f^^iq^) and fsiri^'^) at 
< < Qmax from LCSR. Our main prediction is the integral: 



Hmax Hmax 



AC(0,gLj = :S / d'fpl\fUe)\' = ^tA^ I dq^^^^y^, (2) 



24^3 J - .1 iKblVso J ^ dq 



where Pw = y {f^l + — g2)2/4m^ — is the pion 3-momentum in the i?-meson 
rest frame, and the above equation is valid for £ = e, in the limit mi = 0. As in [12j, 
the value q^ax — 12. GeV^ is adopted. The predicted A(^(0, 12 GeV^) is used to extract 
I Vub I from the most recent B AB AR-collaboration results [21 for the measured partial 
branching fraction integrated over the same g^-region. Furthermore, we predict the 
form factors in the whole semileptonic region by fitting the LCSR results at < q'^ax 
to the z-series parameterization in the form suggested in [l3|. In addition to the q'^- 
distribution of the width, we present a "by-product" observable: the lepton energy 
spectrum in i? — )■ Trii'i. Finally, we calculate the ratio of the B — )• tttUt and B — )• irive 



^Throughout this paper, we assume isospin symmetry for and semileptonic decays and consider 
the B^' -> TT^frui mode for definiteness, denoting it as B — > irEue for brevity. 
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{i = e,fj,) widths, which is independent of \Vub\ and in Standard Model (SM) is fully 
determined by the ratio f^^{q^) / f^^{q^) . We also discuss the current tension between 
the SM prediction and measurements of the leptonic B — )• tv^ width, and suggest to 
use the semileptonic decay B — t- tttUt decay, originating from the same flavour-changing 
interaction, as an additional indicator of new physics. As an illustrative example, we 
consider the influence of a charged Higgs-boson exchange on the B — )• tttv-j- width. 

In what follows, a brief outline of the LCSR method for B ^ ir form factors is given 
in Sect. 2. The choice of the Gegenbauer moments of the pion DA is discussed in Sect. 3. 
In Sect. 4 we present our numerical results for the form factors and for the integrated 
width. In Sect. 5 we discuss the z-series parameterization and present our predictions for 
the whole semileptonic region. In Sect. 6 we discuss the current situation with B — )• rur 
and the decay B — )• tttz/t-, concluding in Sect. 7. 

2 Outline of the method and input 

To obtain the LCSR for the form factors f^^ and /^^ one uses the correlation function 

F^{p,q)=i j <fx e^''-"(7r+(p)|r{n(x)7^6(x),m56(0)i75d(0)}|0) 

= F{q\ {p + qf )p^ + F{q\ {p + qf)q^ , (3) 

of the b ^ u vector current and the i?-meson interpolating current. As explained e.g., 
in [12], the product of the quark operators in the above is expanded near the light-cone, 
provided both external momenta are highly- virtual: {p + q)'^,q'^ *C m'^. The operator- 
product expansion (OPE) result for the invariant amplitudes in ([3|) is obtained in a 
(schematic) form: 

FO^^(<72,(p + g)2) = ^ / VuiT^'HqMp + qf,u„mb,as,fif)^^Hui,fif), (4) 

t 

with a similar expression for F^^^^\ In the above, the pion light-cone DA's (p^^\ui) of 
growing twist t = 2,3,4 are defined as functions of the light-cone momentum fractions 
Ui; for the two-particle DA's ui = u, U2 = 1 — u = u, with the integration over u. 
The DA's are convoluted with the coefficient functions (hard-scattering amplitudes) T*^*-* 
and T^*) at the factorization scale /ij. The currently accessible approximation for the 
light-cone OPE includes the contributions of all two- and three-particle DA's up to the 
twist 4. For the leading twist-2 and twist-3 contributions the coefficient functions are 
calculated in NLO, taking into account the 0{as) gluon radiative corrections. The input 
in the OPE (jH) includes: the 6-quark mass (in the MS scheme) and QCD coupling in 
the coefficient functions, as well as the parameters of the universal pion DA's. 

The dispersion relation for the correlation function ([3]) in the channel of the i?-meson 
interpolating current with momentum p + q is then employed to access the form factors: 

F'^^^iq', {p + qf) = ''^r'ihKl + ••■' 
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fBm% 



m%-{p + q) 



— m 
q2 



1 1 
m-D — m„ 



+fUq')^^^-T^ (6) 



where fs = {B\mJ}ij5d\0) / is the i?-meson decay constant. The ehipses m the 
above relations indicate the integrals over the spectral densities of excited and contin- 
uum i?-states, for which the quark-hadron duality ansatz is used. More specifically, 
one approximates the higher-state contributions in ([5]) and ([6]) by the integrals over 
the spectral density of the calculated invariant amplitudes Im^pj^qyF^^^ {q^ , s) and 

/m(p_)_g)2F'^^^(g^, s), respectively. This approximation brings the effective threshold Sq 
into play. The final form of LCSR is obtained after applying the Borel transformation 
to dS]) and ([6]) replacing the variable (p + q)^ by the Borel parameter M^; e.g., for the 
form factor one obtains: 



= ^ / Im,,^,)^FO^''{q\ s)e-^'^' . (7) 

ml 



The second LCSR obtained from ^ and combined with ([7]) allows one to calculate the 
form factor f%^{q^) ■ Both sum rules are reliable up to q^ax ~ "^b ~ '^^bX^ where x is 
some large scale, independent of mi,, so that at < q^ax truncated light-cone OPE 
can be trusted. 

The interval of M"^ in ([7]) is constrained by combining the two usual criteria for a 
QCD sum rule: smallness of the power corrections (here the contributions of three- 
particle and twist-4 DA's to F^^^), and, simultaneously, a moderate magnitude of the 
hadronic continuum contribution. The interval of Sq is constrained by equating the B- 
meson mass calculated from LCSR to its experimental value. Finally, the decay constant 
fs is calculated from the two-point sum rule with the same Us accuracy. The explicit 
expressions for the amplitudes F'^^^, F^^^ and their spectral functions entering LCSR, 
as well as a detailed description of all pion DA's entering @, can be found in |12|. 

The most important contributions to the LCSR ([7]) originate from the twist-2 and 
twist-3 terms in the OPE (jH). The twist-2 pion DA (/? ^ (ui , ii2 , /x) = fTr(pn{u, fi), is 
normalized to the pion decay constant /tt. The shape of (Pt^{u,^) is determined by 
the coefficients of the Gegenbauer-polynomial expansion (Gegenbauer moments), to be 
discussed in the next section. In the twist-3 pion DA's, the most important input 
parameter is the normalization coefficient = "^^/("T-n + ^^d) related to the quark- 
condensate density. The parameters determining the shapes of the twist-3, 4 DA's are 
known with a sufficient accuracy from the two-point QCD sum rules (see e.g., |15j). 
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3 Gegenbauer moments from the pion e.m. form factor 



For the twist-2 pion DA we use the same approximation as in |12j . 

ipT,{u,Hf) = 6uu(^l + a2(^/)C2^^(ti - u) + al{nf)Cl^^{u - u)j , (8) 

retaining the two nonvanishing Gegenbauer moments and aj and neglecting ah higher 
moments, so that 0^4 = 0. This approximation is justified because the renormalization 
suppresses higher Gegenbauer moments at relatively large scales /if, typical for the LCSR 
([7]). The uncertainties of the input values of 02 4(lGeV) are larger than very small effects 
of the NLO evolution, which we also neglect. 

In [12] these two parameters were constrained by fitting the i? — t- vr form factor 
calculated from LCSR at different to the measured shape of B ^ vr^f^, yielding 
aJ(lGeV) = 0.16 ±0.01 and aJ(lGeV) = 0.04 ±0.01, where the uncertainties only take 
into account the experimental error in the shape. 

Here we refrain from using the B — )• Triiy^ data and obtain an independent constraint 
on the two Gegenbauer moments employing the pion e.m. form factor F^^iQ'^) in the 
spacelike region, defined as 

(vrb ± q)\3r\<p)) = {2p + q),F^{Q^) , (9) 

where fj^ = |u(x)7/,m(x) - \d{x)-i ^,d{x) and = _g2. The LCSR for F^(Q2) derived 
in [TUl [T7] and updated in [TB] is based on the correlation function, similar to Q, with 
virtual u, d quarks instead of the 5-quark and with the axial- vector current instead of the 
pseudoscalar current. This sum rule has NLO accuracy in the leading twist-2 term, with 
the nonleading terms up to twist-6 taken into account. The 0{as/Q'^) term in LCSR 
correctly reproduces the large-Q^ QCD asymptotics of the pion form factor, whereas 
the soft contributions dominate in the intermediate Q'^ region. Importantly, the twist-3 
contribution to the LCSR for F-^(Q^) vanishes in the chiral limit. Altogether, the sum 
rule for the pion e.m. form factor is more sensitive to the twist-2 pion DA than the 
sum rules for heavy-to- light form factors. The LCSR for Ft^{Q'^) with the currently 
achieved accuracy is applicable at intermediate Q^, from O(lGeV^) to a few GeV^; in 
the same region this form factor was accurately measured by the JLab experiment [19j. 
Preliminary comparison of the LCSR with these data was done in [20]. With the pion 
DA given by ([8]) and taking the remaining input from |18j we recalculated the pion 
e.m. form factor at in the region up to a few GeV^ as a function of a2(lGeV^) and 
aJ(lGeV^). As shown in Fig. [TJ the result was fitted to the seven data points at 0.6 
GeV^ < Q2 < 2.45 GeV^, presented in p], yielding: 

aj(l GeV ) = 0.17 ±0.08, aj(l GeV ) = 0.06 ± 0.10 , (10) 

where the uncertainties include experimental errors and the variation of other input 
parameters taken as in |18) . We adopt these intervals for the numerical analysis of the 
sum rule ([7]), neglecting the correlation of the two uncertainties in ()10p . Note that the 
value of 02 presented above is consistent with the direct calculations of this parameter 
in lattice QCD [21] and from QCD sum rules [I5|[22]. 
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Figure 1: The pion e.m. form factor calculated from LCSR as a function of 

Gegenbauer moments a2{l GeV) and aj(l GeV) and fitted (solid) to the ex- 
perimental data points taken from from 



The 7*7 — )• vr transition form factor, where the virtual photon has a spaceUke virtu- 
ahty is another important hadronic matrix element that depends on the properties 
of the pion DA (/^^^(n). The method [23] of combining LCSR with the dispersion rela- 
tion in the photon virtuality predicts this form factor starting from ~ 1 GeV^ . It 
is important that, according to the sum rule approach, the photon-pion transition form 
factor contains a considerable nonperturbative soft contribution. The calculation |23j 
was reconsidered in [20] with the same intervals of 0^,0^ as in [12], revealing a reason- 
able agreement with the data at < 15 GeV^. The most accurate LCSR analysis of 
the 7*7 — )• vr" form factor, including the new twist-6 corrections was carried out recently 
in [H]. It was shown that also the rise of the form factor at high-Q^ observed by the 
BABAR collaboration [2S] (albeit with large errors) can still be accommodated in the 
LCSR prediction at the expense of a moderate "deformation" of the shape of ^Pt^{u). 
Most importantly, one has to increase the coefficient aJ(lGeV) up to ~ 0.2, leaving 
in the ballpark of (fTO]l and (optionally) adding small coefficients Og § to the model. 
We illustrate the numerical influence of such variation of ip-K{u) on the S — )• vr LCSR ([7|) 
in the next section, and confirm the observation made in [23] that the heavy-light form 
factors are only slightly influenced by this modification of Gegenbauer moments. 

The LCSR method can also be used to calculate the D — )• vr form factor. The recent 
analysis in [26] where the same intervals as in [12| were used revealed a very good agree- 
ment with lattice QCD and experiment. We checked that the use of broader intervals 
(jlOp does not produce a noticeable effect, simply because the twist-2 contribution is 
numerically less important in the D — ?■ vr LCSR. 
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4 Numerical results for the form factors and width 



For the numerical analysis of the LCSR ([7]) and the related sum rule for f^^ we slightly 
updated the input used in [12]. First of all, there is practically no change of the b- 
quark mass. According to the last update [27j, we adopt mh{rnb) = 4.16 it 0.03 GeV, 
conservatively inflating the quoted error by a factor of two. As explained in detail in 
[12j . the M5'-mass of the 6-quark is the most suitable mass definition for OPE of the 
correlation function, and the 0{as) contributions to the sum rules are comparably small. 
For the u- and d-quark masses entering the parameter ^.^ = m'^/{mu+md), we follow [2^ 
and use the s-quark mass derived from QCD sum rules |28j and the ChPT light-quark 
mass ratios yielding [mu + md](2GeV) = 8.0 it 1.4 MeV and, correspondingly 

H^{2 GeV) = 2.43 ± 0.42 GeV , (11) 

so that the quark condensate density is {qq){2 GeV) = — (274^^^ MeV)'^. This interval 
is slightly narrower, but remains within the broader range used in |12j. As already 
mentioned in the last section, our intervals for the parameters 05 and aj given in (|10p 
are broader than the ones used in [12]. The rest of the parameters determining the 
nonperturbative objects in the sum rules (DA's and condensate densities), as well as 
the conventions for the renormalization and choice of Os are taken as in [12]. As shown 
there, all nonasymptotic and three-particle contributions of the twist-3 DA's as well as 
the whole twist-4 contribution to LCSR are very small, and the uncertainties in their 
parameters do not produce visible changes in the numerical predictions. 

The "internal" input parameters of our calculation include the renormalization scale 
H, the Borel parameters M and M, the related duality thresholds Sq and Sq in LCSR 
and in the 2-point sum rule for fs, respectively. Here we follow the same strategy as 
in [12], balancing between the smallness of the subdominant contributions (twist-4 and 
twist-2,3 NLO terms) and a reasonable suppression of the integrals over the higher states 
estimated in the quark-hadron duality approximation. The only minor difference with 
respect to the analysis presented in [T2| is that here we stay on a more conservative 
side, allowing for a slightly larger deviation (up to 3%) of the calculated S-meson mass 
from its experimental value. This leads to broader intervals for the Borel parameters 
and duality thresholds. More specifically, we adopt the same default renormalization 
scale /i = 3 GeV as in [12], allowing its variation from 2.5 to 4.5 GeV, and use in LCSR 
the Borel parameter range = (12.0 - 20.0) GeV^ with the threshold parameter 
gradually shrinking from the interval = 37.5 ± 2.5 GeV^ at = 12.0 GeV^ to the 
point Sq = 40.0 GeV^ at = 20.0 GeV^. In the two-point sum rule we vary from 
4.0 GeV^ (s^ = 36.5 ± 2.5 GeV^) to 6.0 GeV^ (s^ = 39.0 GeV^). 

The results of our calculation for f^^{q^) atO<q^< 12.0 GeV^ are shown in Fig. [21 
displaying the separate uncertainties caused by the variation of (a) 02,04, (b) (c) ^u, 
(d) {M'^,Sq} and (e) {M,Sq} within the limits specified above. In addition, in Fig. [2f 
the default central values of the Gegenbauer moments in (jlOp are replaced by a model 
with larger 04(1 GeV) = 0.22 (model III in [21]). The small deviation of the form factor 
remains within our estimated uncertainty due to the 02,04 variation. Since we use a 
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narrow interval of the 6-quark mass from [27], the uncertainties caused by the variation 
of rfib are very smah and not even visible on the plot, hence we do not show them; the 
remaining parameters of DA's and condensate densities generate negligibly small changes 
of the calculated form factors. The sensitivity to the renormalization scale is relatively 
large at approaching g^ax) ^'^^ ^-^so the uncertainties due to the Borel parameter and 
duality threshold are now more pronounced than in [12] due to an enlargement of the 
adopted intervals, whereas the uncertainty due to the twist-3 normalization [i^^ (related 
to the ti, d quark masses and quark condensate) decreases. The numerical results for the 
form factors /bttCq'^) ^'^d /^^(g^) are consistent with what was obtained in [12]; a few 
percent shift of the central value of /^^(O) (presented in Table 1 below) can be traced 
to the modification of the input. 




2 4 6 8 10 12 2 4 6 8 10 12 

tf (GeV^) tf (GeV^) 




2 4 6 8 10 12 2 4 6 8 10 12 

if- (GeV^) (GeV^) 




2 4 6 8 10 12 2 4 6 8 10 12 

if- (GeV^) <^ (GeV^) 



Figure 2: The form factor f^^{q'^) calculated from LCSR with the central input (solid) 
and varying separate input parameters (dashed): (a) aj, (b) ^j^, (c) fi, 
(d) M'^,Sq and (e) M ,Sq; in (f) the result for the model III of Gegenbauer 
moments from is displayed (dashed). 
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Note that the variations of the form factors calculated from LCSR at different are 
strongly correlated, in particular, the shape of the form factor /^^{q^) is correlated with 
its value at q'^ = 0. Quoting separate theoretical errors for each point of accessible q^- 
region, including all these correlations makes the numerical predictions too complex and 
in fact is not necessary, since our main goal is the integrated semileptonic width over this 
region. Instead, we calculate the deviations of this width with respect to the variations 
of individual input parameters, so that the correlations are automatically taken into 
account after the integration over q'^. Our main result is the integral ([2]) calculated 
using the LCSR results for f^^{q^)'- 



AC(0, 12 GeV^ 



4.59 



4,59+1.00 



+0.16 




+0.03 




+0.68 




+0.31 


+0.29 




+0.59 


-0.16 


a2,ai 


-0.03 




-0.46 




-0.39 


-0.47 


A/,so 


-0.32 



















M,so 



ps 



0^85 



(12) 



where the negligibly small uncertainties related to the rest of the input are not shown but 
included in the total error obtained by adding all separate uncertainties in quadrature H. 
Importantly, ()12p has a slightly smaller overall uncertainty than the values of the form 
factor f^^iq"^) at separate g^, due to the abovementioned correlations. 

Using ()12p . we employ the recent BABAR data for the B — )• nivi width. The branch- 
ing fraction integrated from q"^ = to q^ = 12 GeV^ was measured by the BABAR 
collaboration using two different techniques, and the results are: 



Ai3(0, 12 GeV^) = (0.84 ± 0.03 ± 0.04) x lO""^ 
AB{0, 12 GeV^) = (0.88 ± 0.06) x 10"^ 



Taking their weighted average, the total lifetime r^o 
dm) in dg), we obtain: 



2J. (13) 
1.525 lb 0.009 ps and substituting 



IK 



ub\ 



3.50 



-0.38 
-0.33 



±0.11 



th. 



exp. 



X 10 



-3 



(14) 



where the theoretical error corresponds to the estimated total uncertainty in ()12p . 



5 Accessing the large region with z-parameterization 

To extrapolate the calculated form factor, we use the z-series parameterization (see e.g., 
[14^ 130] ) based on the analyticity of the form factors and using the transformation: 



ziq'^,to) 



(niB + m^)2 -q^ - \l (rriB + m^)^ - tp 
{niB + ?7i7r)2 - q2 + y^{mB + m^Y - to' 



(15) 



where to = (w-s + "T-tt)^ — 2y^mBm^y (m^ + m-,^)'^ — q'^^^ is the auxiliary parameter, 
chosen to maximally reduce the interval of z obtained after the mapping (|15p of the 



This replaces our preliminary result A^(0, 12 GeV ) = 4.00lQ g5 quoted in [5] [31 and obtained with 
exactly the same input as in [15], except no data on _B — nlvi were used and broader intervals 
af (1 GeV^) = 0.25 ± 0.15 and aj(l GeV^) = (0.1 ± 0.1) - a^(l GeV^) were adopted. 
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Figure 3: The vector form factor /^^(q^) calculated from LCSR and fitted to the BCL pa- 
rameterization (solid) with uncertainties (dashed), compared with the HPQCD 
(squares) and FNAL/MILC 16] (triangles) results. 



region < q < qinaxi where the LCSR calculation is valid. More specifically, we 
adopt the BCL version [H] of this parameterization, that is, for the vector form factor: 

1 ^ ~ 

fLil') = i_ 2^ E [^(^'' *o)]' • (16) 



As explained in [13], this parameterization has certain advantages with respect to the 
BCL- version |30]. Furthermore, to obey the expected near-threshold behavior, the rela- 
tion 

hN = -^-j^Y.^-lfkK (17) 

fc=o 

is implemented, reducing the number of independent parameters by one. In addition, 
we find it more convenient to keep the form factor at zero momentum transfer /^^(O) 
as one of the fit parameters, correspondingly rescaling the coefficients in the z-series 
expansion. This leads to the same parameterization of the vector form factor as the one 
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Figure 4: The scalar form factor f%^{q^) calculated from LCSR and fitted to the BCL 
parameterization. The notations are the same as in Fig. 0. 



used in [26] 



/l(o) 

1 — jrv? 



B* 



l+Y,bk[z{q^to)''-z{0,toY 



k=l 



-(-1) 



N 



ziq^to)^ -z{0,to) 



N 



(18) 



The scalar form factor f%^{q^) is parameterized in a similar way, except that there is 
no pole factor for the obvious reason: the lowest i?-resonance in the = 0"*" channel is 
located above the B-k threshold. Thus, we use: 



/W) = /L(o) 1 



N 

E 

k=l 



hi 



z{q^,tof -z{Q,to) 




(19) 



where by default 7° JO) = /+^(0). 

We fitted the numerical LCSR prediction for the form factor f^^{q'^) to (jlSp with 
iV = 2 and and fsM^) CS]) with iV = 1, respectively. To increase the "lever arm" 
we also employed the LCSR predictions at negative g^, up to q^i^ ~ ~6.0 GeV^. After 
the mapping (fT5]l . — )• 2; = 0.30 and g^ax = 12 GeV^ — )• 2; = 0.13, so that the 
values of z are sufficiently small to justify truncating the expansion (jl6p . The number 
N of terms in this expansion can be made larger, with no essential change in the fitted 
form factor but with increasing individual uncertainties for the coefficients bk and with 
large correlations between them. We checked that the upper bounds on the expansion 
coefficients bk following from the OPE of the two-point correlation function of the vector 
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(GeV') 



Figure 5: (colour online) The normalized q'^ -distribution in B ^ ttIv obtained from 
LCSR and extrapolated with the z- series parameterization (central input- solid, 
uncertainties -dashed). The experimental data points are from BABAR: (red) 
squares 12], (blue) triangles /5/ and Belle f^J: (magenta) full circles. 



bj^u currents (see [H] for detailed expressions) are far from being saturated for low N. 

We also used the analogous bounds for the scalar form factor obtained recently in [31] • 

Altogether, the OPE bounds play a role starting from N = 5. 

The fitted values of /^^(O) = /^7r(0) and of the slope parameters 6i, ft? are presented 

in Table [H together with the numerically important uncertainties, the latter revealing 

significant correlations. With these results we extrapolate the form factors at > q'^ax 

and compare with the lattice QCD results. This is shown in Fig. [3] for f^^iq"^) and in 
Fig.Hfor 

The dashed curves in these figures are obtained by adding separate variations of the 
form factors in quadrature at each g^, so that the variations of the solid curves corre- 
sponding to the central input are bounded within the area between the upper and lower 
dashed curves. As expected, the uncertainties of the form factors extrapolated to larger 
g^, exceed the ones calculated at smaller q^. This circumstance, however, does not play 
a significant role for the integrated widths, since the integration over the phase space 
suppresses the semileptonic width in the large g^-region. 

Our predictions for f^^ are, within errors, in a reasonable agreement with the lattice 
QCD results obtained by HPQCD [5] and Fermilab/MILC [6j collaborations. We also 
observe an agreement with the normalization and shape of the form factors obtained by 
the QCDSF collaboration 0, in particular, they predict /^^(O) = 0.27 ± 0.07 ± 0.05. 
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Parameter 


centr. value 


{a2,a4} 








{m\so} 




0.281 


+0.002 
-0.003 


+0.018 
-0.014 


-0.005 
+0.008 


+0.010 
-0.022 


-0.010 
+0.016 


^1 


-1.62 


+0.43 
-0.44 


-0.06 
+0.05 


+0.53 
-0.07 


+0.30 
-0.49 




6? 


-3.98 


+0.56 
-0.57 


-0.28 
+0.23 


+0.96 
-0.08 


+0.28 
-0.42 





Table 1: Fitted parameters for z-series parameterization of the form factors f^^ {q^) 
and their uncertainties due to the variations of the input parameters. 

Furthermore, we estimate the total width of — )• irivn in units of l/|Kif,P and the 
integral ([2]) for the large g^-region: 

— ^r(S ^ 7:^ui) = AC(0,26.4 GeV') = 7.71^^-^^ ps-^ , 

\Vub\ 

AC(16 GeV^ 26.4 GeV^) = 1.88l[]|^ ps"^ . (20) 

Our prediction for the latter integral has to be compared with the lattice QCD results 
presented below, in Table [H 

In Fig. [5] we plot the predicted (7^-shape in — t- ttIv^ obtained by calculating the 
normalized differential width {l/T)dT / dq^ from LCSR at < < q^ax from the 
z-series parameterization at g^ax < 9^ ^ {f^B — 'n-Tr)'^- The estimated uncertainties are 
naturally smaller than in FigO because the variations of the form factor normalization 
cancel in the ratio of the differential and total widths. Our result is compared with the 
measured (/^-distributions. We use the partial Ai3 spectrum obtained by the BABAR 
collaboration in 6 bins and 12 bins from the two independent analyses [5] and |3j, re- 
spectively. For the normalization we employ the corresponding central values of the 
measured total branching fractions: B{B^ -f vr-^+z^^) = (1.41 ± 0.05 ± 0.07) x 10"^ [2] 
and B{B^ vr^^+i/^) = (1.42±0.05±0.07) x 10"^ [3j. The analogous 13-bin distribution 
measured by the Belle collaboration [4J and normalized by their total branching fraction 
B{B^ TT~l+V() = (1.49 ± 0.04 ± 0.07) x 10-"^ is also shown. Fig. [5] reveals a general 
agreement of our prediction for the g^-shape and experimental results, however, only 
within still large uncertainties of both experiment and theory. In particular, the shape 
of the form factor fitted from the BABAR data in pj to the same BCL parameterization 
with two parameters yields: (6i/6o)babaj? = — 0.67 ± 0.18, whereas our result for the 
same ratio is (bi/bQ)LCSR = -l-10to:27 • 

A "byproduct" observable that was not yet measured in i? — )• vr^f^, is the distribution 
of the lepton energy in the S-meson rest frame shown in Fig. [6] in the normalized form. 
In the electron or muon semileptonic i?-decay, at a given lepton energy, the distribution 
dT{B — )• ■Klvi)/dEi contains the integral of \fg^{q^) p over the region < (7^ < ImsEi 
(the expression for this distribution, also at 7^ 0, can be found, e.g., in [32j). Hence 
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the correlations between the normahzation and shape of the form factor somewhat reduce 
the uncertainties in this distribution. It also has a more pronounced slope than the 
distribution. 




0.0 0.5 1.0 1.5 2.0 2.5 

E, (GeV) 



Figure 6: Lepton energy spectra for B — vr^f^ at m£ < Eg < (m^ + mj) / {2mB) for 
i = fi,T. Solid (dashed) lines correspond to the form factors calculated at the 
central input (indicate the uncertainties). 



6 5 — > ti/t and B — > ttti/t 

Currently, the leptonic width B — )• measured by both BABAR and Belle collabora- 
tions (see Table 2) is larger than the SM prediction: 

BiB- ^ Ti?r) = ^\Vub\'m',mB (l - ^) flrj,- , (21) 
Svr V 

if one employs fs predicted from lattice QCD or QCD sum rules, together with \Vub\ 
extracted from B — )• vrff^. The recent discussions on this situation are mostly con- 
centrated on the value of \Vub\- Indeed, the tension decreases, if one uses in ([2T]) the 
somewhat larger value of \Vub\ extracted from the inclusive b ^ u decays. On the other 
hand, the CKM fits \33\ |3l] yield a smaller \Vub\, consistent with the determinations 
from B — )• Tvii^i. 

Let us emphasize that, independent of the actual \Vub\ value, there exists a tension 
between the ratio of semileptonic and leptonic B widths and the QCD predictions for 
the two relevant hadronic matrix elements /^^((7^) and /b- To demonstrate that, we 
define the following observable: 

o 2. ^ ^BB^^iuMh^i) ( iB_\ ^ Ac(gf,gi) 

^s//Wi,g2J- B{B^TVr) \tb^) (Gy8vr)m2mB(l-m2/m|)2/2 ' 
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where the partial branching fraction AB and the integral defined as in ([2]), are taken 
over the same region q\ < tf' < of the momentum transfer. 

The above equation for the ratio Rgji follows solely from the V — A structure of the 
weak currents in SM and Vub cancels out in the ratio. The form factor fg^ and decay 
constant /s entering r.h.s. are obtained by one and the same QCD method: lattice 
QCD or the combination of LCSR and QCD sum rule. In Tables [2] and [3] we collect the 
inputs for this equation, obtained from different measurements and QCD calculations. 
The disagreement between the calculated and measured ratio Rgji goes beyond the 
theoretical and experimental errors, especially in the case of the lattice calculations 
which have smaller uncertainties. 



Exp. 


A5(10-^) [Ref.] 


B{B ^ TUr){l^~^) [Ref.] 




BABAR 
Belle 


0.32 ± 0.03 [2| 
0.33 ±0.03 ±0.03 [3] 

0.398 ± 0.03 [2 


1.76 ±0.49 [371 [38] 


0.281°:^? 


QCD 


AC(ps-i) [Ref.] 


/B(MeV) [Ref.] 


Rs/i 


HPQCD 
FNAL/MILC 


2.02 ±0.55 i5j 


190 ± 13 [35] 
212 ±9 [36] 


0.52 ±0.16 
0.46 ±0.10 



Table 2: The ratio Rgji for the region 16 GeV^ < < 26 A GeV^, measured and calcu- 
lated from ()22p using the lattice QCD results. The weighted average over the 
two BABAR measurements is taken and all errors are added in quadrature. 

Decreasing further the theoretical and experimental errors in (j22p . especially in the 
B — )• Ti/r width, becomes therefore a very important task. Possible effects beyond the SM 
in B ^ TVr are already being discussed in the literature, and, in particular, B Dtu^ 
is proposed as a channel which has common new physics contributions with the leptonic 
B decay (see e.g., [ID] and references therein). 

Here we would like to attract attention to another semileptonic channel: B — )• tttUt, 
although it is experimentally very demanding. Earlier this channel was discussed e.g., 
in |32[ 141] . Note that this channel has the same combination of quark and lepton 
flavours as 5 — )• tv^- In the SM, the B — ?■ tttz^t- decay differs only kinematically 
from the semileptonic modes with the muon or electron. A convenient, Vub-independent 
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Exp. 


AB{10-*) [Ref.] 


B{B TUr){W-^) [Ref.] 


Rs/i 


BABAR 


0.88 ±0.06 [2] 
0.84 ± 0.03 ± 0.04 [3] 


1.76 ±0.49 [371138] 


0.52+0-?o 

— U.lz 


QCD 


AC [Ref.] 


/B(MeV) [Ref.] 


Rs/i 


LCSR/QCDSR 


4- 59^0:85 [this work] 


210 ± 19 02] 


n 07+0-28 



Table 3: The same as in Table [2] for the region 

< g2 < 12.0 GeV^ where the QCD sum 

rule results are used. 



observable Il3] is the ratio 



dT{B —?- -KTUrMdq^ _ {q^ - m-r)^ ( _|_ 



dr{B ^ TTii^e)/dq^ (g2)2 V 2g2 

3m2(m|-m2)2 |/0jg2 )| ^ 
4(m2+2g2)m|p2 |y+j^2 



where i = e or fi and is neglected. It is determined by the ratio of the scalar and 
vector i? — 7- vr form factors, hence, it has a somewhat smaller uncertainty than the 
individual form factors. The ratio (j23|) is plotted in Fig. [S] in the kinematically allowed 
region < q"^ < {niB — JTItt)^, using our predictions for the form factors. It grows at 
large q'^, due to the kinematical suppression of the vector channel contribution at small 
Ptj. The strong correlations between the vector and scalar form factors calculated from 
LCSR with one and the same input result in a small uncertainty in their ratio. The 
tau-lepton energy spectrum for B — )• ttti' shown in Fig. [6] is another observable which 
depends on both vector and scalar form factors. 

To investigate the influence of new physics on both leptonic and semileptonic B decays 
with a r-lepton, we include in the effective Hamiltonian of the b — t- utVt- transition an 
intermediate charged Higgs-boson contribution adopted in the same generic form as in 

m ■■ 



Gf ( 



- "^T^ ^ [ds +9P75jbT{l - 75)j^r| ± h.c. , (24) 

The admixture of new physics in i? — )• ti/j- and B — )• ttti/t is then determined, respec- 
tively, by the pseudoscalar and scalar parts of the new interaction. In particular, the 
leptonic width ()2ip gets multiplied by (1 — gp)"^; and therefore, with this choice, the 
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3.5 
3.0 
2.5 
2.0 
1.5 
1.0 
0.5 
0.0 



TTT V. 



dr (B -* 71 fi Vp) / dq 




(GeV') 



Figure 7: (colour online) Ratio of differential decay widths, defined in \23\) (solid), with 
shaded (green) area indicating the uncertainties. Also shown is the effect of 
adding a charged Higgs-hoson contribution with gs = —0.4 (dashed, red) and 
gs = 2A(dash-dotted,red) 



B —7- TVr width vanishes at gp = 1. Accordingly, the r.h.s. of (|22p acquires a factor 
1/(1 — gp)"^- Also the ratio ()23p is modified by multiplying the scalar form factor with 
an additional factor: 

fUl')-^(l-'-^)fBM% (25) 

The addition of the new interaction can fill the gap between the calculated and ex- 
perimentally measured ratio i?^/; in ()22p if one allows for gp ^ 0. Taking, e.g., the sum 
rule prediction for this ratio from Table 3, and adding the new physics contribution, 
R^H — )■ Rg/i/{l — gp)'^, we equate it to the experimental value and find that gp ^ 
is allowed within one of the following two intervals: gp = —(0.4 it 0.2\th ± 0.2\exp) or 
gp = 2.4 lb 0.2|j/j lb 0.2|exp. Assuming the parameter gs in the same ballpark as gp (in 
fact, gs = gp in MSSM, however, only with positive values), we display in Fig. [7] the 
modified ratio (|23p . adding the new physics contribution with gs = —0.4 and gs = 2.4. 
In fact, the second option may already be excluded by the B — DtUt analysis (see e.g., 
[S]). In addition, in Fig. [H the ratio of total semileptonic widths is plotted as a function 
of gs- Note that the deviation due to new physics can be larger than the uncertainty 
due to the hadronic form factors. 
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Figure 8: Ratio of total semileptonic decay widths into r and fi(ore) (solid), with shaded 
(green) area indicating the uncertainties shown as a function of the parameter 
gs determining the charged Higgs-boson coupling in B ^ iirVr 



7 Discussion 

The precise determination of the CKM matrix element Vub is mandatory for stringent 
tests of the quark-flavour content of SM. Due to recent work in lattice QCD and in QCD 
sum rules, combined with constraints from analyticity and unitarity, the values for \Vi^b\ 
extracted from B — t- T^iv are becoming quite precise. The main result of this paper is the 
prediction for the partial width of this decay in the g -region < g2 < 12 GeV^ in terms 
of which is expressed as a weig hted integral AC(0, 12 GeV^) over the squared 

form factor f^^{(f') ■ This integral has smaller uncertainties than the values of the form 
factor at separate and allows the \ Vub\ extraction with an accuracy approaching 10%, 
somewhat better than in the previous LCSR analysis [12] where only the value of the 
form factor at = was used. We extract a value for \Vub\ using recent BABAR data 
[21 [3]. We hope that also the Belle collaboration will in future provide the integrated 
width in this region. 

We also employed the z-series parameterization with BCL-ansatz and extrapolated the 
LCSR form factors to the whole kinematic region. The B ^ it form factors f^^{q^) and 
f^^{q'^) obtained from LCSR are, within comparable uncertainties, in agreement with 
the recent lattice QCD results and with the measured g^-shape of the B — )• Trtui width. 
A combined fit of the form factors calculated at small q'^ from LCSR and at large from 
lattice QCD to a common z-series parameterization is possible (see e.g. the previous 
analysis in [H]) but is beyond the scope of this paper. 

A further improvement would need more precise measurements of the shape to con- 
trol the input parameters used in LCSR. On the theoretical side, the renormalization 
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scale dependence can be reduced by including NNLO corrections to the hard scatter- 
ing amplitudes, and also by separating the renormalization and the factorization scales. 
Improving the duality approximation is more difficult and demands a knowledge of ra- 
dially excited states in B channel. Another perspective is a simultaneous global fit of 
three different LCSR's to the data on i? — t- Triui and on the pion electromagnetic and 
transition form factors, with scanning over the allowed region of input. 

The value for \Vub\ obtained from B — )• nii^i is somewhat lower than what is extracted 
from inclusive B decays, as well as the one extracted from B — )• rP, but the significance 
is still too small to be conclusive. On the other hand, our result is completely compatible 
with the results from the CKM fits ^331 ED- 

The impact of the recent measurements oi B ^ tv on Vub has not improved the 
situation concerning our knowledge of this quantity. However, the ratio of leptonic and 
semileptonic widths is independent of Vub and may either be regarded as a test for a 
possible non-standard contribution (like, e,g. a charged Higgs-boson exchange) or as 
a test of our understanding of QCD. In turn, a lack of understanding of /s and the 
form factor fg^ will severely limit the sensitivity to a "new physics" contribution. It is 
interesting to note that both QCD sum rules as well as lattice calculations tend to yield 
larger values of a suitably defined ratio of semileptonic versus leptonic widths. This 
means that we have either a problem in our understanding of QCD matrix elements or 
that there is really a substantial new-physics contribution in i? — )• rP, which also makes 
the channel B — )• iirvr very interesting. A significant test of these statements has to 
await more data, in particular on B ^ tD. 
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